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ABSTRACT 


The  thermomechanical  response  of  a  uniaxial  bar  with  thermovis coplastic 
constitution  is  predicted  herein  using  the  finite  element  method.  After  a 
brief  review  of  the  governing  field  equations,  variational  principles  are  con¬ 
structed  for  the  one  dimensional  conservation  of  momentum  and  energy  equations. 
These  equations  are  coupled  in  that  the  temperature  field  affects  the  displace¬ 
ments  and  vice  versa. 

Due  to  the  differing  physical  nature  of  the  temperature  and  displacements, 
first  order  and  second  order  elements  are  utilized  for  these  variables,  respec¬ 
tively.  The  resulting  semi-discretized  equations  are  then  discretized  in  time 
using  finite  differencing.  This  is  accomplished  by  Euler's  method,  which  is 
utilized  due  to  the  stiff  nature  of  the  constitutive  equations. 

The  model  is  utilized  in  conjunction  with  stress-strain  relations  devel¬ 
oped  by  Bodner  and  Partom  to  predict  the  axial  temperature  field  in  a  bar  sub¬ 
jected  to  cyclic  mechanical  end  displacements  and  temperature  boundary  condi¬ 
tions.  It  is  found  that  spacial  and  time  variation  of  the  temperature  field 
is  significantly  affected  by  the  boundary  conditions. 

TABLE  OF  SYMBOLS 


t  -  time 

P  -  axial  internal  resultant  force 
p  -  axial  externally  applied  force  per  unit  length 
x  -  axial  coordinate  dimension 
a  -  axial  stress  component 
A  -  cross-sectional  area 

T  -  end  traction  in  units  of  force  per  unit  area 
x 

s  -  surface  area 
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Table  of  Symbols  (cont.) 


S  -  area  of  the  longitudinal  surface  of  the  bar 
c 

e  -  axial  strain  component 
u  -  axial  displacement  component 

a  -  internal  state  variable  representing  axial  inelastic  strain 
E  -  Young's  modulus  in  the  axial  coordinate  direction 

a  -  coefficient  of  thermal  expansion  in  the  axial  coordinate  direction 
T  -  temperature 

T  -  reference  temperature  at  which  no  deformation  is  observed  at  zero  load 
R 

-  internal  state  variable  representing  drag  stress 

q  -  axial  component  of  heat  flux 

k  -  coefficient  of  axial  thermal  conductivity 

C  -  specific  heat  at  constant  elastic  strain 
v 

p  -  mass  density 

r  -  internal  heat  source  per  unit  mass 
L  -  length  of  the  bar 

INTRODUCTION 

It  is  well  known  that  mechanical  and  thermodynamic  coupling  are  signif¬ 
icant  in  metallic  solids  [1-11].  The  author  has  recently  developed  a  model 
capable  of  predicting  this  coupling  effect  in  thermovis coplastic  metals  [12]. 
In  the  previous  paper  a  cyclic  strain  control  loading  on  a  sample  of  INI 00  at 
1005 °K  (1350°F)  was  used  to  predict  a  temperature  rise  of  approximately  3.7°K 
per  cycle  when  the  strain  amplitude  was  2%  and  the  specimen  was  adiabatically 
insulated. 

The  focus  of  the  current  research  is  to  consider  the  effect  of  thermal 
boundary  conditions  on  this  same  process.  The  introduction  of  these 
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conditions  causes  the  strain  and  temperature  fields  to  be  inhomogeneous  even 
though  the  stress  field  is  homogeneous  if  the  bar  is  prismatic.  This  spacial 
variation  in  the  field  variables  causes  the  process  to  be  difficult  to  model 
because  the  thermomechanical  constitutive  equations  are  highly  nonlinear  stiff 
differential  equations.  In  this  paper  the  finite  element  method  is  utilized 
*  to  spatially  discretize  the  dependent  variables  displacement  and  temperature, 

and  the  finite  difference  method  is  employed  for  timewise  discretization. 

This  process  results  in  a  set  of  highly  nonlinear  algebraic  equations. 

Since  the  thrust  of  this  research  is  to  obtain  accurate  results  without 
regard  to  numerical  efficiency,  the  results  are  obtained  via  the  relatively 
inefficient  but  accurate  method  of  simply  utilizing  successively  smaller  time 
steps  along  with  refined  spatial  mesh  to  obtain  a  convergent  and  therefore 
accurate  solution  for  the  temperature  and  displacement  fields  both  spatially 
and  as  a  function  of  time  for  a  cyclically  imposed  end  displacement. 

The  physical  interest  in  the  problem  is  to  determine  the  effect  of 
temperature  boundary  conditions  on  the  predicted  temperature  rise  in  a  bar  sub¬ 
jected  to  cyclic  mechanical  loading.  It  is  found  from  the  analysis  that  the 
introduction  of  these  nonadiabatic  boundary  conditions  causes  significant  axial 
temperature  gradients.  Since  nonadiabatic  conditions  cannot  be  avoided  in 
experimental  research,  it  is  concluded  that  experimental  tests  of  this  type 
should  be  viewed  with  caution  when  their  purpose  is  to  construct  constitutive 
relations. 

PROBLEM  SOLUTION 

Field  Problem  Description 
The  following  field  equations  are  given: 
a)  equilibrium  [13], 

3P  =  -p  (x)  ,  O) 

9x  X 
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where  the  axial  resultant  P  is  defined  by 


c 


b)  strain-displacement  relation 


c)  thermomechanical  constitution, 

a  =  E[e  -  aj  -  a(T  -  XR)]  .  , 

(£ ,  T ,  , a j )  ,  i  -  1  »z  , 


and  (2) 

(3) 

(4) 

(5) 

and  (6) 


q  * 


(7) 


where  z  is  the  total  number  of  internal  state  variables;  and 

d)  conservation  of  energy 

[<Ee  -  EOl  +  EaTR)  +  EO2  1  ||]  -  Ea  1  ||  -  pcv  ||  -  |§  +  pr  -  0  .  (8) 

The  conservation  of  mass  is  satisfied  trivially  and  the  second  law  of  thermo- 
dynamics  has  been  previously  shown  to  be  satisfied  by  the  above  equations  [14- 
lb].  It  should  be  noted  that  equilibrium  equation  (1)  satisfies  equilibrium 
in  the  axial  coordinate  direction  only  in  an  average  sense  over  the  cross-section. 
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The  above  6+Z  equations  (excluding  definition  (3))  define  a  nonlinear 
initial-boundary  value  problem  (together  with  appropriate  thermal  and  mech¬ 
anical  initial  and  boundary  conditions)  in  which  the  following  dependent  vari¬ 
ables  are  sbught  as  functions  of  x  and  t:  0,  £,  u,  q,  T,  P,  and  a  . 

For  convenience  the  domain  is  defined  to  be  of  length  L,  so  that  boundary 
and  initial  conditions  are  of  the  form: 


u  (x ,  0)  =  Uq  =  known 
T  (x ,  0)  =  Tq  =  known 


initial  conditions 


(9) 


and 


essential 

boundary 

conditions 


/  u(0,t) 
u(L,t) 
|  T  (0, t) 
^  T(L,t) 


=  u^  =  known  or  P(o.t) 
=  u^  =  known  or  P(L,t) 
=  =  known  or  q(0,t) 
=  =  known  or  q(L,t) 


=  P^  =  known 
L 

=  P  =  known 
=  q^  =  known 
=  qk  =  known 


natural 
boundary 
conditions.  (10) 


It  is  now  assumed  that  a  =  a(x)  so  that  equation  (2)  reduces  to 


P  =  aA 


(ID 


Therefore,  substituting  (4)  into  (5)  and  this  result  into  (11)  gives 


P  =  EA 


[# '  “i  -  “<T  -  v] 


(12) 


The  above  result  is  now  substituted  into  (1)  to  obtain 


3_ 

3x 


EA 


[If  -  al  -  0(1  -  Tr)]  1  -  - 


px(x) 


(13) 
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which  represents  the  differential  equation  relating  displacements  and  temper 
ature  to  the  applied  load  p  (x) . 

Equations  (4)  and  (7)  are  next  substituted  into  energy  balance  law  (8) 
and  this  result  is  integrated  over  the  cross-sectional  area  A  to  obtain 


where  it  has  been  assumed  that  all  field  variables  depend  on  x  and  t  only. 
Now  define 


Q  = 


dA  =  -k 


3T 

3x 


A 


(15) 


Careful  inspection  of  equations  (13)  and  (14)  will  indicate  that  these 
equations ,  together  with  internal  state  variable  growth  laws  (6)  and  initial 
and  boundary  conditions  (9)  and  (10),  represent  a  well-posed  boundary  value  prob¬ 
lem  in  terms  of  the  2+z  dependent  variables  u,  T,  and  <x  . 

Solution  Procedure 

The  field  problem  is  to  be  solved  analytically  using  the  semi-discretized 
finite  element  technique  with  timewise  finite  differencing.  In  order  to  ac¬ 
complish  this,  differential  equations  (13)  and  (14)  must  first  be  written  in 
a  suitable  variational  form. 

Variational  Equations 

Consider  first  equation  (13).  This  governing  equation  is  integrated  against 
a  suitably  smooth  test  function  v  =  v(x)  over  the  domain  of  some  element  : 


6 


*e+i  r 


Ot,  - 


a(T  -  T  ) 
R 


■] 


+  p 


X 


dx  =  0 


(16) 


Integrating  by  parts  results  in 


re+l  r  1 

■  J  ai "  a(T '  vjdx  = 


v  EA 


[s  -  “i  -  “«  -  v] 


e+1 


x 


X 


X 

/ 


e+1 


v  p  dx 
rx 


(17) 


Substituting  equation  (12)  into  the  boundary  term  thus  results  in 


-  f  “•  If  [if  -  “i  -  a<1  -  v]  dx  ■ 

X  ~ 


"v(3Cen)  P(xe+1)  +  v(xe)  P(xe) 


3^  e+1 

X  ~ 


v  p  dx 
x 


(18) 


Now  consider  equation  (14).  Once  again  the  governing  equation  is  inte¬ 
grated  against  a  suitably  smooth  test  function  w  =  w(x)  over  the  domain  of  the 


element  ft  : 

e 


re+1  f  r 

J  “  r  l^1 


|e|H-  Eaj  +  EaIR)  ^  +  Ea2  T  || 


-AEcff  jt3x  -  ApCy  St  +  A  jx 


(«)■ 


+  A  pr 


r  dx  =  0 


(19) 
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Integrating  the  heat  flux  term  by  parts  results  in 


(20) 


where  equation  (15)  has  been  substituted  into  the  boundary  terms. 

Finite  Element  Spacial  Discretization 

Quadratic  displacement  and  linear  temperature  fields  are  now  chosen  within 
each  element: 


3 

u(x,t)  =  2  ue  ipB  ,  x  <  x  <  x  .  ,  and  (21) 

.  ,  i  i  e  e+1 

1=1 


T(x,t)  =  2  T®  <|>J 

i=l 


X  <  X  <  X  , 

e  e+1 


where  u®  =  u®(t)  and  T®  =  T®(t)  are  the  nodal  displacements  and  temperatures, 
respectively,  and  =  ^?(x)  an<^  $7  =  4h(x)  are  quadratic  and  linear  shape 
functions,  respectively  [17]. 

Appropriately,  v  and  w  are  endowed  with  the  properties  of  u  and  T,  respec¬ 
tively,  so  that 
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(23) 


v  =  ipf  i  =  1,2,3 


w  5  <f>?  i  -  1,2 


Substitution  of  equations  (12)  and  (21)  through  (23)  into  variational  principle 
(18)  results  in 


“I  I-  -  Ol*  ■  '  P<W  +  *i(xe)  P<Xe) 


j.e+1 

-  /  *: 

■V 


p  dx  i  *  1,2,3 

X 


(24) 


The  above  may  be  written  in  the  following  compact  form 


I  Ke.  ue  +  Z  S6.  Te  =  Fe  ,  i  =  1,2,3, 

3=1  1J  J  3=1  J  1 


(25) 


where 


x 


K6. 
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e+1  j  i  e  j  |  e 

d^.  dip 

EA  —  -t— ^  dx 
dx  dx 


i  =  1,2,3;  j  =  1.2,3  ; 


(26) 
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x 

/ 


e+1  ,  ,e 

dip. 

EAa  j— —  dx 

dx  Tj 


i  =  1,2,3;  j  =  1,2 


(27) 
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X 


6+1  d*f 


/dt|)7 

“dT  (-“l 

X 


+  aTt>)dx 
K 


x  i 

/e+1 

X 


Px<ix  , 


i  =  1,2,3. 


(28) 


Similarly,  substitution  of  equations  (21)  through  (23)  into  equation  (20)  re¬ 
sults  in 


x 


/ 

X 


e+1 


-“arkl  £  TM)  +  A*i 


(*  ^  (Ji 


+  EaTR  +  Ea 


-  Ea 


(j,  i 


e ,  e 

4> 

mm 


dx  = 


X 

e+1 

(2 

\ 

+  / 

y 

\ 

e,  e  \ 

T  <j>  J  -  pr 
mm  / 

X 

e 

m 

/ 

dx. 


i  =  1,2, 

Equations  (29)  may  be  written  in  the  following  form: 


(29) 


3  2  e+1 

I  K6.  ue  +  E  S&.  Te  +  f  (pf  A 

.  ,  ii  l  .  ,  it  i  J  i 

3=1  J  J  3=1  J  x 


Ea‘ 


I  Te  <pe 
J-l  3  3 


2  dT 
v  m  ,  e 
2  TT  4> 


m=l 


dt  rm , 
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-E 


+  pr 


dx 


dTe 
_ m 

dt 


m 


e+1 


3a, 


-  "  /  A(-Ea1  +  EaTR)  -  Q(x.) 

x 


i  -  1,2, 


(30) 


where 


*e+l  dip .  3a, 

/  “  »i  s1  ir dx 

x 

e 


i  =  1,2;  j  =  1,2,3;  and  (31) 


*e+l  d<J>?  d<j>e 

S6.  =  -  /  kA  ”j  —  -r-1  dx  ;  i  =  1 , 2;  j  =  1 , 2.  (32) 

ii  J  dx  dx 

j  x 

e 


Finite  Difference  Timewise  Discretization 

Time  dependence  in  equations  (6)  and  (30)  is  handled  via  finite  differ¬ 
encing,  Although  higher  order  approximations  may  be  used,  Euler  forward  dif- 

e  e 

ference  approximations  are  now  entered  for  the  time  rate  of  change  of  a^,  T^, 

and  ue. 
m 

3a? 

(x ,  t )  s  [aR  (x,t  +  At)  -  a£  (x,  t)  ]  /At,  k  =  1,  ...,z  (33) 

dTe 

3-2  (t)  s  [Te  (t  +  At)  -  T  (t) ] /At,  m  =  1,2  ,  and  (34) 

at  m  m 

due 

-5-2  (t)  s  [ue  (t  +  At)  -  u  (t)]/At,  m  =  1,2,3-  (35) 

ctt  m  m 


ii 


Substitution  of  (33)  through  (35)  into  finite  element  equations  (30)  gives 


3  2 

Z  K?.  ue  +  Z  S?.  T 
j=l  13  3  j-1  13  3 


e+1 


/ 


Ea 


X 


Z 

j-1  3  3, 


u 

m=l  x 


T  (t  +  At) 
m 


At 


<(t)V 

m 


■) 


-Ea 


j-1  J  J 


3  /u*(t  +  At)  -  u*(ty,|£ 


m-1' 


At 


9 


3x 


-PC 


2  /T  (t  +  At) 
m 


if 

m=l  ' 


Tm(t) 


At 


) 


'm 


+  Pr 


>  dx 


x  .  i 

/e+1 

A 4»J  [-Ea1 


3a.. 

(t)  +  EaTR]  (t) 


-Q<xi> 


i  =  1,2 


(36) 


The  above  may  be  written  as  follows: 


Z  K? .  uf  +  I  S6.  Te 

j-i  J  j*i  J 


2  2  2 
+  Z  Z  C.  ..  Te  T.e  +  Z  D.  .  T. 
k=l  j-1  ljk  3  k  j-1  1J  3 


23  2 
+  Z  Z  E.,  .  T?  Ue  +  Z  G.,  t: 
k=l  j-1  lk3  k  3  j-1  13  3 


+  Z  H. .  T.  =  FT 

J-1  13  3  1 


i  =  l,2. 


(37) 
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where 


te+1 

:ijk  =  /  A*lzf  ♦*  dx 


i  =  1,2;  j  =  1,2;  k  =  1,2  ,  (38) 


x 


-r 


ti'ux^)2  +  dx. 


i  -  1,2  ,  (39) 


X 


Di2  E 


■•/ 


e+1 


A<j>?  —  [Tj(t)4>j4>j  +  T®(t)(<|>*)2]  dx. 


i  -  1,2  >  (40) 


Em  .  = 
ikj 


x  .  i 

•/ 


_  3^e 

.  ,e  Ea  ,e  rj  , 

A+i  IF  *k  IT  dx 


i  =  1,2,;  k  =  1,2;  j  =  1,2,3  ,  (41) 


X 


Gil 


■/ 


e+1 


A<|> 


e  Ea 
i  At 


drbf  dip*  dip* 

u*(t)  <Dj^  +  «;ct)  ^  ^4-u3e(t)  ^ 


dx,  i  =  1,2  ,  (42) 


x 


’i2 


/*e+^  „  r  3n»;  1 

E  /  Ul(t)*2  dT  +  u2(t)Cp2  dT  +  U3(t)<p2  3F-Jdx  *  i  = 


x 


1,2  ,  (43) 


H.  .  = 
13 


■/ 


e+1 


,e  pCv  . e  , 

A<h  dx 


i  =  1,2;  j  =  1,2; 


x 


and  (44) 


Fe  = 

i 


?e+l 

-  /  < 
X  " 


pC  /  2 


At 


2  T?(t)<f>?  /+  pr 

j  =  l  J  3' 


dx 


/e+1 

M )®  [~Ea1 

A 


3a1 

(t)  +  EaTR]  (t)  dx 


-Q(Xi)  i  =  1,2 


(45) 
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Equations  (37)  may  be  written  equivalently  as  follows: 


u6  +  Z  SB.  Te  =  Fe 
J  J-l  1J  J  1 


where  K^f .  and  F®  are  as  defined  previously,  and 


2 

K® .  +  Z 
k=l 


2 

Z 

k-1 


C...  T.  + 
vjk  k 


'ij 


+  G.  . 


+  H.  . 


(46) 


and  (47) 


(48) 


The  above  equations  may  be  adjoined  with  equations  (25)  to  obtain  the  following 
set  of  nonlinear  equations  for  each  element. 


5x1 


(49) 


where  all  nonlinearity  is  contained 


in  [S],  {F®},  and  {F6}. 


Global  Assembly  and  Boundary  Conditions 

Global  assembly  is  accomplished  in  the  standard  way  using  the  Boolean 
matrix  [17].  Interelement  continuity  is  guaranteed  by  setting 


+  pf 1  -  0 

9 

and  (50) 

+  Ilif'1  -  0 

• 

(51) 

Boundary  conditions  are  implemented  in  the  standard  way: 
boundary  conditions  are  handled  by  placing  one  on  the  diagonal 


1)  essential 
of  the 
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appropriate  row  and  zeros  off  diagonal  in  the  stiffness  matrix,  and  the  speci¬ 
fied  value  of  the  essential  variable  on  the  right  hand  side;  and  2)  natural 
boundary  conditions  are  implemented  directly  to  the  right  hand  side. 

Solution  of  the  Nonlinear  Algebraic  System 

Initial  conditions  are  used  for  the  first  time  step.  The  time  step  At 
is  supplied  for  each  load  increment  and  boundary  conditions  are  incremented 
directly  from  supplied  input  functions. 

The  internal  state  variable  is  handled  in  equations  (23)  and  (45)  by 
using  equations  (35).  is  initialized  according  to  reference  18.  The  non¬ 

linear  stiffness  matrix  [S]  is  initialized  using  nodal  temperatures  and  displace¬ 
ments  from  the  previous  time  step.  The  displacements  and  temperatures  at  time 
t+At  are  then  estimated  directly  and  without  iteration  by  utilizing  equations  (49) 
for  very  small  time  steps. 
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EXAMPLE  PROBLEMS 


In  order  to  completely  define  an  example  problem  it  is  necessary 
to  specify  internal  state  variable  growth  laws  (6).  Numerous  models 
have  been  proposed  for  crystalline  metals  [18,19].  Since  it  is  not  the 
purpose  of  this  research  to  compare  these  models,  a  relatively 
established  model  proposed  by  Bodner  and  Partom  [20]  has  been  chosen. 
This  model  contains  two  internal  state  variables:  the  inelastic  strain 
(a  )  and  the  drag  stress  (c^).  The  growth  laws  for  these  variables  are 

given  by 


a  i 


]71  exp 


(52) 


and 


a  2 


m(Zi  -  a2)  aai  -  AZi 


(53) 


where  DQ>  n,  m,  Z  ,  Z  ,  and  r  are  experimentally  determined  material 
constants . 

For  the  purpose  of  modeling  the  temperature  gradient  in  a  specific 
component,  a  hypothetical  problem  has  been  chosen  using  material 
properties  representative  of  Inconel  100  at  1005°k  (1350°F).  The 
material  and  geometric  properties  are  given  in  Table  I.  The  geometry  is 
representative  of  a  cylindrical  uniaxial  bar  which  is  2,50  inches  long 
and  0.25  inches  in  diameter.  It  has  previously  been  shown  that  Bodner 
and  Partom1 s  model  accurately  predicts  the  stress-strain  behavior  of 
INI  00  under  uniaxial  loading  conditions  for  both  monotonic  and  cyclic 
strain  controlled  loadings  [12,18]. 
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pCv  =  5.032  MPa/°K 

a  =  13.14  x  10  ^  in/in/°K 

k  =  21.0  x  10-6  MPa  M2 

sec°K 

E  =  146.86  x  103  MPa 
A  =  7.12557  x  10'5  m2 
T_  =  1005 °K 

K 

L  =  .06350  m 

3 

Dq  =  10  x  10  in/in 

n  =  0.70 
m.  =  2.57 
Zx  =  1015.0 

Zj  =  600.0 

r  =  2.66 


Table  I-»  Material  Properties  for  IN  100  at  1005°K  (1350°F) 
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Utilizing  the  material  properties  described  above,  the  following 
effects  have  been  studied  using  the  model  developed  herein: 

1 )  the  effect  of  variation  of  strain  rate  on  the  time  dependence 
of  temperature  at  the  midpoint  of  a  monotonically  extended  uniaxial  bar 
which  is  insulated  on  the  longitudinal  boundaries  [Figs.  1-4]; 

2)  the  spacial  variation  of  temperature  for  the  case  above  [Fig. 

5]; 

and 

3)  the  effect  of  end  temperature  boundary  conditions  on  the 
temperature  at. the  center  of  a  uniaxial  bar  which  is  held  at  fixed 
temperature  at  the  end  points  and  subjected  to  cyclically  imposed  end 
displacements  [Figs.  8  and  9]. 

Examples  1  and  2  are  constructed  primarily  to  determine  the  effects 
of  thermomechanical  heating  on  the  stress-strain  behavior  of  uniaxial 
constitutive  specimens.  It  is  found  in  examples  1  and  2  that  if  a 
specimen  is  mounted  in  an  experimental  apparatus  which  has  massive  grips 
simulating  a  fixed  temperature  boundary  condition  there  can  be 
substantial  axial  temperature  gradients  induced  in  a  time  dependent 
boundary  layer  near  the  ends  of  the  specimen.  On  the  other  hand,  these 
boundary  conditions  do  not  appear  to  substantially  affect  the  predicted 
stress-strain  behavior,  especially  when  the  strain  measurement  is  taken 
between  the  thermal  boundary  layers  near  the  grips.  Therefore,  it  would 
appear  that  the  standard  procedure  for  obtaining  stresses  and  strains  in 
uniaxial  bars  is  not  substantially  affected  by  mechanically  induced, 
axial  temperature  gradients  when  the  grips  are  at  fixed  temperature 
equivalent  to  the  initial  specimen  temperature  and  the  bar  is  loaded 
monotonicalaly.  However,  it  should  be  noted  that  massive  grips  which 
are  mounted  outside  a  furnace  could,  by  their  much  lower  temperature 
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STRESS  (MPA) 


1400.0 


Fig.  1. 


Predicted  Stress  vs.  Strain  for  a  Uniaxial  Bar 
Pulled  at  Various  Constant  Strain  Rates. 


TENSION 


ABSOLUTE  STRAIN  (x-L/2) 


Predicted  Temperature  vs.  Absolute  Strain  for  Monotonic 
Deformation  Histories  Described  in  Fig.l.  (e=  +.0142  sec 


TEMPERATURE  CHANGE  (°K) 


ABSOLUTE  STRAIN  (x-L/2) 


Fig.  3.  Predicted  Temperature  vs.  Absolute  Strain  for  Monotonic 

Deformation  Histories  Described  in  Fig.  1.  (e=  +.00142  sec-  ) 


2 


TEMPERATURE  CHANGE  (°K) 


Predicted  Temperature  vs.  Absolute  Strain  for  Monotonic 
Deformation  Histories  Described  in  Fig.  1.  (e=  +.000142  sec  1) 


than  the  initial  speciman  temperature,  may  induce  significant  error  in  predicted 
strains  if  the  strain  is  measured  by  dividing  relative  displacement  by  some 
gage  length. 

The  final  example  demonstrates  that  under  cyclic  loading  conditions  the 
above  conclusions  may  not  necessarily  be  true,  especially  when  the  specimen  is 
subjected  to  high-cycle  fatigue  and  at  high  strain  rates.  There  is  definitely 
a  trend  towards  an  increasing  mean  temperature  in  the  bar,  and  this  mean  tempera¬ 
ture  is  strongly  affected  by  the  thermal  boundary  conditions  as  well  as  the 
loading  rate.  Although  it  would  be  interesting  to  determine  the  mean  temperature 
rise  in  a  cyclic  fatigue  test,  the  current  algorithm  precludes  this  analysis  due  to 
the  extremely  large  computer  times  necessary  to  predict  only  a  few  cycles  of 
response  (approximately  43.8  GPU  minutes  on  an  Amdahl  470/V6  for  the  example 
demonstrated  in  Figs.  6  and  7). 

Example  3  also  demonstrates  another  interesting  phenomenon  which  may  be 
significant  in  large  space  structures.  If  the  bar  is  perfectly  insulated  the 
mean  temperature  rise  per  cycle  for  the  relatively  slow  loading  rate  shown  in 
Fig.  6  is  3 . 7°K,  whereas  if  the  ends  of  the  bar  are  held  at  a  fixed  temperature  of 
1005°K,  the  mean  rise  is  1.0°K  per  cycle.  Faster  loading  rates  show  less  difference 
between  the  adiabatic  result  and  the  fixed  end  temperature  result.  Since  many 
of  these  structures  are  expected  to  be  extremely  flexible  truss— like  configurations, 
a  typical  metallic  member  which  undergoes  some  yielding  (which  might  be  desirable 
in  order  to  induce  natural  damping)  might  in  fact  undergo  substantial  enough 
heating  during  vibrational  response  such  that  the  material  properties  could  be 
adversely  affected,  thus  resulting  in  a  material  related  failure  of  the  structure. 
However,  further  investigation  is  needed  on  this  last  issue  since  it  is  expected 
that  the  primary  form  of  heat  flux  off  of  space  structures  will  be  via  radiation 
on  the  longitudinal  surfaces  of  the  truss  member.  Since  the  current  analysis 
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has  treated  these  surfaces  as  insulated,  no  general  statements  can  be  made  at 
this  time  regarding  thermomechanical  heating  in  space  structures. 


24 


TEMPERATURE  CHANGE  (°K) 


STRESS  (MPa) 


Stress-Strain  Curve  and  Strain  input  curve  for  cyclic  load 


STRAIN  <cw/cm> 


CONCLUSION 


The  current  research  has  attemped  to  demonstrate  the  effects  of 
mechanical  loading  on  one-dimensional  temperature  gradients  in  a  class 
of  viscoplastic  media.  Due  to  the  nonlinearity  and  stiffness  of  the 
field  equations,  it  was  necessary  to  utilize  a  numerical  algorithm. 
This  algorithm  has  been  shown  to  be  very  inefficient  for  solving  even 
one-dimensional  examples.  Therefore,  it  is  apparent  that  significant 
refinement  of  the  procedure  will  be  necessary  before  multi-dimensional 
analyses  can  be  performed  by  this  method.  Specifically,  it  would  be 
significant  to  determine  the  effect  of  transverse  temperature  gradients 
on  the  stress-strain  behavior  of  constitutive  specimens.  Furthermore, 
the  effects  of  thermal  boundary  conditions  on  the  longitudinal  surface 
needs  attention.  The  author  is  currently  studying  a  perturbation 
technique  for  more  efficient  solution  of  these  issues. 

The  above  points  notwithstanding,  the  current  research  demonstrates 
some  important  results.  These  are: 

1 )  The  axial  temperature  gradient  in  a  viscoplastic  uniaxial  bar 
is  strongly  affected  by  the  thermal  boundary  conditions  on  the  ends. 

2)  The  end  temperature  boundary  conditions  can  cause  temperature 
gradients  which  are  substantial  enough  to  induce  spacial  variations  in 
stress  and  strains  which  invalidate  the  standard  procedure  of  using 
average  quantities,  although  when  grips  are  mounted  within  a  furnace  at 
spacially  constant  temperature,  it  appears  that  the  standard  procedure 
is  accurate. 
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3)  There  is  a  trend  toward  increasing  average  temperature  in 
cyclically  loaded  bars;  whether  or  not  this  effect  is  significant  is 
strongly  dependent  on  the  thermal  boundary  conditions  and  the  loading 
rate. 
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